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Using a first-principles parameterized lattice-gas Hamiltonian we study the adsorbate ordering 
behavior at atomic steps of a Pd(lOO) surface exposed to an oxygen environment. We identify a 
wide range of gas-phase conditions comprising near atmospheric pressures and elevated temperatures 
around 900 K, in which the step is decorated by a characteristic O zigzag arrangement. For catalytic 
processes like the high-temperature combustion of methane that operate under these conditions our 
approach thus provides first insight into the structure and composition at a prominent defect on the 
working surface. 



INTRODUCTION 

As prominent defects at solid surfaces atomic steps are 
commonly perceived as playing some kind of special, if 
not decisive role for the surface properties or function in 
materials science applications. When it comes to the in- 
teraction with a reactive environment, steps are viewed 
as providing particularly active sites which can have a 
pronounced effect on the catalytic activity or act as nucle- 
ation centers for adsorbate-induced morphological tran- 
sitions like oxide formation or corrosion in our oxygen- 
rich atmosphere. When aiming to qualify this role at 
the atomic scale an important first task is to identify the 
structure and composition at the step edge under real- 
istic gas-phase conditions. On the modeling side, the 
first-principles atomistic thermodynamics approach [l[ 
has brought considerable progress to this end for ideal 
low-index surfaces. In its prevalent form, this approach 
compares the stability of a variety of structural models 
in contact with a given gas-phase reservoir. Extending it 
to the study of steps depends then only on the ability to 
perform the underlying electronic structure calculations 
for corresponding structural models. For metal surfaces, 
a viable route to this is to resort to supercell geometries 
for vicinal surfaces, i.e. to compute periodic arrays of 
steps that can be suitably cast into a periodic boundary 
condition framework. Unfortunately, the computational 
cost connected to such inherently large supercell calcula- 
tions severely limits the accessible surface unit-cell sizes 
and total number of structural models that can be com- 
puted. This prevents to date a proper exploration of the 
huge configurational space of possible step structures and 
allows at best for approximate treatments of entropic ef- 
fects at elevated temperatures. 

In the present study we address these limitations by 
suitably combining the thermodynamic approach with 
a first-principles parameterized lattice-gas Hamiltonian 
(LGH) [l|, [3, 13| • The concept is illustrated with the ap- 
plication to the on-surface oxygen adsorption at Pd(lOO), 
and there specifically to the decoration of close-packed 
(111) steps in an oxygen atmosphere. We employ large- 



scale density-functional theory (DFT) calculations for 
corresponding low-index and vicinal surfaces to deter- 
mine the lateral interactions between adsorbates at the 
(100) terrace sites and at sites close to or at a (111) 
step. In a first stage the resulting multi-site LGH is 
then used to generate a large pool of possible structural 
models with varying on-surface O content. Within the 
thermodynamic framework the stability of these mod- 
els is evaluated for a wide range of gas-phase conditions 
and compared to known more complex 0-induced recon- 
structions at the low-index surface, namely an ultrathin 
(VS X \/5)R27° (henceforth simply termed \/5) surface 
oxide and thick bulk- like PdO films 4, 5, 6]. Intrigu- 
ingly, this identifies a distinct range of O chemical po- 
tentials where a characteristic zigzag step decoration is 
most stable. Explicitly accounting for the configurational 
entropy we confirm with grand-canonical Monte Carlo 
(MC) simulations that this ordered structure is indeed 
stable up to quite elevated temperatures, in fact up to 
much higher temperatures than the concurrent ordered 
adsorbate structure at the Pd(lOO) terrace. At near at- 
mospheric pressures such temperatures around 900 K are 
representative for catalytic applications like the high- 
temperature combustion of methane 0], for which our 
study thus provides first insight into the structure and 
composition at a prominent defect on the working sur- 
face that may serve as a basis for ensuing kinetic studies, 
e.g. using kinetic Monte Carlo simulations ^li]. 

THEORY 

The atomic-level energetics for our first-principles sta- 
tistical mechanics approach is obtained from DFT cal- 
culations, which we perform within the full-potential 
augmented plane wave (LAPW/APW+lo) framework Q 
and using the generalized gradient approximation for the 
exchange-correlation interaction In preceding work 
exactly the same computational setup [lO| was validated 
and used to systematically investigate the on-surface ad- 
sorption of oxygen at Pd(lOO) and the series of 
Pd(lliV) vicinals (with = 3, 5, 7), which exhibit (111) 
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FIG. 1: Schematic view of a Pd(lOO) surface with (111) steps. Shown are the five distinct adsorption sites considered in the 
first-principles lattice-gas Hamiltonian: The fourfold coordinated (100) hollow sites at the terraces (The), at the upper (Thu) 
and lower (Thl) rim of the step, as well as the two threefold (111) hollow sites at the step edge, one with a coordination to two 
step atoms (Sh2) and one with a coordination to one step atom (Shi). Additionally represented are four characteristic ordered 
structures discussed in the text: A p(2 x 2) and a c(2 x 2) adsorbate phase at the (100) terrace, a zigzag structure where the 
step is decorated with O atoms in alternating Thu and Sh2 sites, and a trilayer structure where all Thu and Sh2 are occupied 
(Pd = large spheres, Pd step atoms = large dark spheres, O = small spheres). 



steps and (100) terraces of increasing width p^]. In the 
limit of low O coverages the energetically most favorable 
adsorption sites obtained were the fourfold coordinated 
hollow sites at the extended Pd(lOO) surface (The), as 
well as the two threefold coordinated hollow sites (Shi 
and Sh2) offered by the (111) step. Furthermore, only 
the adsorption energetics into the (100) hollow sites di- 
rectly at the upper and lower rim of the (111) step (Thu 
and Thl) was noticeably different from that at the corre- 
sponding The site at the extended Pd(lOO) surface. We 
correspondingly focus our present study on these five 
distinct sites as depicted in Fig. [H and describe the 
site-specific adsorption of oxygen atoms on the stepped 
Pd(lOO) surface with a lattice-gas model, in which any 
system state is defined by the occupation of these sites in 
a lattice. The total free binding energy of any configura- 
tion (comprising static and vibrational contributions) is 
expanded into a sum of discrete lateral interactions be- 
tween adsorbates at the lattice sites. We consider pair 
interactions between any combination of sites up to dis- 
tances that correspond to the third nearest neighor (NN) 
shell at Pd(lOO), as well as the most compact many- 
body interactions in form of trio interactions between 
any three site combination within a second NN shell dis- 
tance. Judging from our preceding systematic analysis 
of lateral interactions at the extended Pd(lOO) surface 

ll\ we are confident that this truncated LGH covers the 
ordering behavior at the accuracy level of interest to this 
study, and we will comment on this accuracy level and 
steps to further validate the truncation below. 

In order to generate a LGH with predictive accuracy 
the unknown static on-site and lateral interaction ener- 
gies are parameterized with the DFT binding energies 
computed for a set of ordered configurations, closely fol- 
lowing the recipe detailed for the low-index surface in ref. 

ll|. In brief, we determine those on-site energies and 
interactions influenced by the (111) step from large su- 
percell calculations for the Pd(117) vicinal, and those of 
undisturbed terrace sites from calculations for the low- 



index Pd(lOO) surface. In all cases, the ordered struc- 
tures were fully geometry optimized, which as already 
detailed in ref. [12l | leads to pronounced relaxations of up 
to 20% of the first interlayer spacing at the step ridge at 
high adsorbate coverages. The four on-site energies for 
the Shi, Sh2, Thu and Thl sites are accordingly taken 
from calculations with one O atom inside a Pd(117)- 
(1 X 3) surface unit-cell, and the on-site energy for the 
terrace The site from calculations with one O atom in- 
side a Pd(100)-(3 X 3) surface unit-cell. With distances 
to the periodic adatom images that are in both cases al- 
ways larger than 8.37 A (corresponding to the sixth NN 
shell at Pd(lOO)) we verified that these calculations give a 
very good approximation to the low-coverage limit. The 
remaining total of 70 lateral interaction energies are de- 
termined by least-squares fitting to DFT binding energies 
of 102 ordered configurations described in similarly sized 
supercell geometries. 

A first indication for the reliability of the resulting 
LGH expansion comes from the computed leave-one-out 
cross validation score H , [13, 14 1 of 29 meV / adatom, 



which provides a measure for the average accuracy with 
which the truncated LGH can predict the energetics of 
adatom configurations not included in the fitting proce- 
dure. The consistency of the expansion procedure is fur- 
thermore reflected by the similarity of equivalent LGH 
parameters that can be independently determined by fit- 
ting to the low-index or the vicinal DFT data. In all cases 
such interaction energies differed by less than 20meV 
and we verified that choosing either or the other value 
did not affect any of the conclusions presented below. 
Due to the large mass mismatch between O and Pd, the 
vibrational contribution to the LGH parameters can be 
approximately considered through its leading term which 
arises from the change of the O2 gas-phase stretch fre- 
quency to the 0-substrate stretch mode Exten- 



sive test calculations for various ordered configurations 
involving all five adsorption sites revealed only a small 
coverage-dependence of this term, which is correspond- 
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ingly modeled as constant shifts of the five on-site ener- 
gies by the zero point energy changes computed in the 
low-coverage limit. 

RESULTS 

In a first stage we exploit the low computational cost 
to evaluate the thus established first-principles multi-site 
LGH to extensively explore the huge configuration space 
and identify low-energy ordered step structures. This 
is achieved by simulated annealing MC runs in periodic 
boundary condition simulation cells containing (60 x 40) 
(100) surface unit-cells and one (Iff) step along the 
shorter cell axis. For fixed and initially random O cover- 
ages ranging up to a half filling of all terrace hollow sites 
10® MC steps are used to continuously quench from an 
initial temperature of 2000 K. Depending on the O con- 
tent, these simulations unanimously lead to patches of 
increasing size of four distinct ordered structures that are 
schematically shown in Fig. [TJ At the lowest coverages 
the O adatoms decorate alternating Thu and Sh2 sites 
at the (III) step in a characteristic zigzag structure. At 
increasing coverages domains of the experimentally well 
characterized p(2 x 2) overlayer [3| start to form at the 
(100) sites next to the upper edge of the fully zigzag dec- 
orated step. Only once the O content in the simulation 
cell exceeds the one required to fully cover the entire ter- 
race with the p{2 x 2) structure, does a further decoration 
of Thu and Sh2 step sites occur. However, the formation 
of the resulting periodic step structure which we hence- 
forth term trilayer, cf. Fig. [T] goes then directly hand in 
hand with the growth of a c(2 x 2) arrangement [3] at the 
terrace sites. The first-principles LGH thus correctly re- 
produces both experimentally characterized ordered ter- 
race structures. As already demonstrated in ref. [lH this 
holds also for other aspects of the low-coverage phase di- 
agram of the ideal Pd(IOO) surface hke the order-disorder 
transition, which is better described than by any previ- 
ous empirical LGH. Without existing experimental data 
such empirical modeling would furthermore currently be 
precluded for the adsorbate ensemble near the (III) step, 
whereas the present first-principles approach clearly pre- 
dicts two new ordered structures in which the O adatoms 
decorate the step. 

For the identified structural models we now account 
for the effect of a finite gas-phase environment through 
the atomistic thermodynamics approach Here, the 
surface is considered to be in equilibrium with an oxy- 
gen gas-phase reservoir characterized by a chemical po- 
tential A/io(T, j>) summarizing the two-dimensional de- 
pendence on pressure p and temperature T and using 
molecular O2 as reference. In its prevalent form this ap- 
proach neglects the effect of configurational entropy at 
the surface and then simplifies to approximately com- 
puting the Gibbs free energy of adsorption for ordered 
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FIG. 2: Computed Gibbs free energy of adsorption AG^'^ for 
O adsorbate structures at a (111) step at Pd(lOO) and using 
the clean surface as zero reference. Using the nomenclature 
defined in Fig. [l]the labels refer to the structure at the step 
and at the neighboring terrace, e.g. zigzag-|~p(2 x 2) stands for 
a zigzag step decoration coexisting with the p(2 x 2) adsorbate 
ordering at the terrace. Additionally indicated by differently 
shaded backgrounds are the stability regions obtained for the 
ideal Pd(lOO) surface in previous work [13, [Si, namely those 
of the p(2 X 2) overlayer, of the \/5 surface oxide structure 
and of a thick bulk-like PdO film 



structural models at any given A/ io( T, p) from their cor- 
responding binding free energies |I5l . Il6l | . This allows to 
readily compare the stability of most distinct structural 
models over a wide range of gas-phase conditions, and is 
thus ideally suited to also assess under which gas-phase 
conditions the surface structure and composition at the 
step is appropriately described by our on-surface LGH. 
We correspondingly employ this approximate approach 
here not only to obtain a first understanding of the sta- 
bility of the presented step models, but also to contrast 
this with the stability ranges of more complex 0-rich 
phases identified previously for the ideal Pd(IOO) surface 
[E H IHj IHI • Specifically, these are the V5 surface 
oxide structure, which corresponds to a sub-nanometer 
thin film of PdO(IOI) on the surface and bulk PdO 
to represent thick bulk-like oxide films. The results are 
summarized in Fig. [5]and reveal a large stability range of 
the zigzag decoration of the (III) step, which over most 
of this range coexists with the p{2 x 2) adsorbate phase 
on the (100) terrace. The trilayer arrangement becomes 
only more favorable, i.e. exhibits a lower Aff*"^, at quite 
elevated chemical potentials and coexists then directly 
with the denser c(2 x 2) overlayer on the terrace. In light 
of the also shown stability ranges of the oxidic phases 
at ideal Pd(IOO) we can therefore identify quite a range 
of gas-phase conditions (— I.43eV < A^o < — I.I6eV) 
that are not yet 0-rich enough to induce oxide forma- 
tion and where the zigzag step decoration is not just due 
to kinetics, but instead represents a thermodynamically 
stable phase in its own right. In contrast, the elevated 
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FIG. 3: Order parameters 'i' sensitive to the order of the 
zigzag step and p(2 x 2) terrace structure, as determined 
by grand-canonical MC simulations for p = 10"'^ atm. At 
this pressure, the temperature range shown corresponds to 
— 1.43eV < Afio < — 1.16eV, where the approximate ther- 
modynamic results of Fig. [2]predict the stability of the zigzag 
structure or the coexistence of zigzag + p{2 x 2) as indicated 
by the differently shaded backgrounds. Equivalent results are 
obtained for simulations at p = 10~^ and 1 atm, with varia- 
tions of the critical temperatures of the order of 100 K. 



ably employed to determine the critical temperature Tc 
for the order-disorder transition at fixed oxygen pressure 
Interestingly, the disturbance exerted by the step 
on the ordering behavior on the terrace is found to be 
rather short ranged on the accuracy level of interest to 
this study. Already four atomic rows away from the im- 
mediate step edge the computed critical temperatures are 
to within ±50 K identical to those obtained at the ideal 
Pd(lOO) surface. Figure [3] therefore only contrasts the or- 
dering behavior at the immediate step edge with that at 
ideal terrace sites for p = 10~^ atm. Whereas the approx- 
imate thermodynamic approach predicts the presence of 
the p{2 X 2) ordered structure at the terrace up to a tem- 
perature of 980 K (corresponding to A/ito = — 1.37eV at 
this pressure), the low order parameter obtained in the 
grand-canonical MC simulations shows that the terrace 
adlayer has in fact long disordered at these temperatures. 
This underscores that a full understanding of the surface 
structure and composition at such high temperatures can 
only be obtained through the explicit evaluation of the 
partition function as enabled by the first-principles LGH 
approach. In contrast, the Tc = 950 K deduced from the 
inflection point of the order parameter for the zigzag dec- 
orated step validates that this ordered structure prevails 
indeed up to surprisingly high temperatures. 



A/iQ {T, p) at which the trilayer arrangement at the step 
becomes more favorable than the zigzag decoration fall 
well into the stability range of the oxidic surface recon- 
structions. This suggests an interpretation of the trilayer 
structure in form of a kinetic precursor in the oxide for- 
mation process, and in this respect it is intriguing to note 
that its geometric structure coincides in fact with one row 
of the \/5 surface oxide 5] . 

At near atmospheric pressures the identified stability 
range of the zigzag step decoration (— 1.43eV < A^o < 
— 1.16eV) corresponds to rather elevated temperatures 
of the order of 900 K [l^. At such temperatures, con- 
flgurational entropic effects can not be uncritically dis- 
missed. Considering the high melting temperature of Pd 
and the high kink formation energy of the close-packed 
(111) step, we expect these effects to be dominated by a 
possible disorder in the adsorbate ensemble. We there- 
fore scrutinize the insight provided by the approximate 
thermodynamic approach by fully accounting for such ef- 
fects with grand-canonical MC simulations based on the 
first-principles LGH. In these simulations, the effect of 
the (111) step on the adsorbate ensemble in its vicin- 
ity is evaluated by separately monitoring the order and 
coverage for each row of sites parallel to the step edge 
using the same simulation cells as for the simulated an- 
nealing runs. Since the prevalent ordered structures in 
the targeted (T, p)-range are the p{2 x 2) at the terrace 
and the zigzag decoration at the step, order parameters 
that are sensitive to these lateral periodicities can be suit- 



SUMMARY 

Through a combination of first-principles statistical 
mechanics techniques we have been able to identify a wide 
range of oxygen environments, in which (111) steps at 
Pd(lOO) will be decorated by a characteristic zigzag oxy- 
gen structure. This includes near atmospheric pressures 
and temperatures around 900 K, i.e. gas-phase conditions 
that are representative for an important catalytic process 
like the high-temperature combustion of methane. Due 
to the low energetic cost of the close-packed steps, this 
defect will be a frequent structural motif at the surface 
with corresponding potential influence on the function 
in such applications. In this respect it is interesting to 
notice that the same ordered zigzag O structure was re- 
cently characterized at the (111) steps of Pd(lll) vicinals 
in the same range of o xyg en chemical potentials, albeit 
at lower temperatures [21|. Using our LAPW/APW-l-lo 
DFT setup we compute initial and final state O Is core 
level shifts [l^ that are identical to within 40 meV for the 
O atoms adsorbed at the geometrically equivalent Sh2 
step sites in the two structures. Since core level shifts 
are a sensitive probe of the local bonding properties, this 
suggests highly similar reactivities of the oxygen atoms 
adsorbed at the decorated (111) steps at both prominent 
facets of Pd nanoparticles. If these oxygen atoms play 
a determining role for the catalytic function, this simi- 
larity would provide an intriguing atomic-scale interpre- 
tation for the enigmatic structure insensitivity reported 
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recently for the methane combustion at these elevated 
temperatures Q- While our study focuses on establish- 
ing a first methodological access to determine the struc- 
ture and composition of a prominent defect like atomic 
steps at surfaces exposed to realistic environments, this 
immediately exemplifies the far-reaching insight that this 
kind of first-principles statistical mechanics approach can 
provide. 
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